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Abstract 

We introduce a two-dimensional lattice model for the description of knot- 
ted polymer rings. A polymer configuration is modeled by a closed polygon 
drawn on the square diagonal lattice, with possible crossings describing pairs 
of strands of polymer passing on top of each other. Each polygon configura- 
tion can be viewed as the two-dimensional projection of a particular knot. We 
study numerically the statistics of large polygons with a fixed knot type, us- 
ing a generalization of the BFACF algorithm for self-avoiding walks. This new 
algorithm incorporates both the displacement of crossings and the three types 
of Reidemeister transformations preserving the knot topology. Its ergodicity 
within a fixed knot type is not proven here rigorously but strong arguments 
in favor of this ergodicity are given together with a tentative sketch of proof. 
Assuming this ergodicity, we obtain numerically the following results for the 
statistics of knotted polygons: In the limit of a low crossing fugacity, we find 
a localization along the polygon of all the primary factors forming the knot. 
Increasing the crossing fugacity gives rise to a transition from a self-avoiding 
walk to a branched polymer behavior. 
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The statistical properties of self-avoiding polymers are now very well understood, 
mainly thanks to the famous equivalence of the problem to a field theory, leading to 
a natural derivation of scaling laws and critical exponents In three dimensions 
however, closed polymers will generally form knots, which, for self-avoiding objects, 
will restrict the accessible configurations. In the last few years, some progress has 
been made towards answering some basic questions about knots, showing for instance 
that a sufficiently long polymer is knotted with probability one [Q, [3], or computing 
the distribution of random knots [[| . However, the statistics of a closed polymer with 
a fixed knot-type is less well understood. Indeed, the usual field theory approach 
does not account for this distinction of the knot type, but corresponds rather to a 
summation over all the possible knot topologies. In particular, the natural question 
of how the critical exponents depend on the knot topology remains an open issue. 
Knotted configurations are encountered in the description of closed DNA molecules, 
with an apparent influence of the knot type on some of their properties J5J . 

An attempt to understand the role of a fixed knot topology was done in || [7|, || 
|S], TO [TT[] where a numerical study of knotted lattice polygons in the cubic lattice Z 3 
was performed. There it was shown that, while the connectivity constant and the size 
exponent v for the radius of gyration are independent of the knot type, the entropic 
exponent a for the number of accessible configurations clearly depends on the knot 
at hand. More precisely, a seems to depend only on the number of primary knots 
(factors) and increases by one for each added factor in the knot factorization. 

In practice, all the characterizations of knots involve only their projection in a 
two-dimensional (2D) plane. Therefore, nothing prevents us from considering a two- 
dimensional model for knotted polymers. Such a model should be in principle much 
simpler to study, both analytically and numerically. The aim of this paper is precisely 
to introduce and study numerically a particular model of "projected" polygons on a 
2D lattice, describing 2D knotted polymers. In a 2D model for knotted polymers, 
the self-avoiding constraint in two dimensions is released to allow for what we call 
"crossings" . By crossings, we actually mean two strands of the chain passing on top 
of each other in the projection. The object remains self-avoiding in the sense that 
the underlying and overlying strands cannot be exchanged by passing through each 
other. With this restriction, each closed polymer has a fixed knot type, which can be 
preserved in the deformations. Still the object is two-dimensional in the sense that its 
entropy corresponds to 2D deformations only. In particular, we will recover in some 
limits some of the usual exponents of 2D self-avoiding walks (SAW). 

The paper is organized as follows. In the first section, we introduce our lattice 
model of projected self-avoiding polygons. Section 2 describes the local elementary 
moves which we use to deform the polygon and explore its accessible phase space 
at a fixed knot topology. The issue of ergodicity is also discussed in this section 
where arguments for a proof of ergodicity within a fixed knot type (but not a full 
proof) are given. The elementary moves are then performed according to a Monte 
Carlo Metropolis algorithm, drawn to reproduce the most natural grand canonical 
statistical weight, with a fugacity K per bond of polygon, and a fugacity w per 
crossing. The updating rules are discussed in section 3. Several improvements of the 
algorithm are described in section 4. They involve both non local deformations and 
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Figure 1: An example of a PSAP configuration on the square diagonal lattice with 
six crossings and the topology of a trefoil knot. 



multiple chains run in parallel. Our numerical results are presented in sections 5, 6 
and 7 for w — 1, w — > and w > 1 respectively. We gather our conclusions in section 
8. 

1 The Model 

We present here the model of Projected Self-Avoiding Polygons (PSAP's) that we 
shall use to describe two-dimensional knotted polymers. The model is defined on the 
square diagonal lattice in two dimensions, namely the simple two dimensional square 
lattice Z 2 completed with the diagonals of the squares. We treat each diagonal as 
an elementary edge of the lattice, i.e. we consider that the intersection point of the 
two diagonals is not a vertex of the lattice. The bonds of the polymer can sit on 
all the edges of the lattice, either vertical, horizontal or diagonal, with at most one 
bond per lattice edge. Different rules apply for the edges and vertices of the simple 
square lattice on the one hand and for the diagonal edges on the other hand. We 
impose a strict self-avoidance at the vertices of the square lattice, i.e. we do not allow 
two parts of the polymer to either cross each other or even to touch each other at 
these vertices. We allow crossings to take place inside a square at the intersection 
point of its two diagonal edges. More precisely, we use the diagonal edges only for 
crossings, i.e. we impose that a diagonal edge can be occupied if and only if the 
perpendicular diagonal edge in the same square is also occupied, with the two bonds 
on these edges forming a crossing. For each pair of occupied diagonal edges inside a 
square, we distinguish between two different possible crossings according to which of 
the two bonds of polymer passes on top of the other. We can thus view crossings as 
a two dimensional projection of a pair of bonds in three dimensions, with one bond 
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Figure 2: The three types (I, II and III) of Reidemeister moves on the projections of 
knots. Reversing all the crossings also corresponds to allowed transformations. 



lying on top of the other. Finally, we limit ourselves to closed polymers. Figure [I] 
shows an example of an allowed configuration with six crossings and the topology of 
a trefoil knot. 

Our model can be seen as a simple extension of the usual model describing self- 
avoiding polymer loops on the square lattice, also refered to as Self- Avoiding Polygons 
(SAP). The new ingredient here is the possibility for the polymer, which is strictly 
self-avoiding on the square lattice, to have crossings taking place on pairs of perpen- 
dicular diagonal edges inside a square. Again, these "crossings" are viewed as the 
two dimensional projection of two bonds of polymer passing on top of each other, 
with thus two distinct allowed configurations according to which of the two bonds is 
on top of the other. Thanks to this distinction, our model still describes a particular 
self-avoiding object in the sense that it mimics the projection in two dimensions of a 
polymer which would be self-avoiding in three dimensions. We will refer to our model 
as a model of Projected Self- Avoiding Polygons (PSAP's). Since we use closed poly- 
gons, a PSAP will in general form knots, i.e. will be the two-dimensional projection 
of a three dimensional knot. In this sense, we can speak of our model as describing 
two-dimensional knotted polymers. 

To fully specify the model and to study the corresponding statistics, we need to 
assign to each PSAP configuration its weight. Since our aim is the study of knotted 
polymers of a fixed knot type, we want to attach a non-zero weight only to those 
configurations which have the desired knot-topology. We will be mainly interested 
below in rather simple topologies, i.e. that of the unknot (0), that of the trefoil 
knot (3i), that of the figure eight knot (4 X ) and that of the composite knots made 
of two trefoils (3i#3i) or a trefoil and a figure eight (3i#4i). In practice, we will 
start with an initial configuration fixing the knot type and we will explore the phase 
space accessible by performing successive transformations (moves) on the PSAP which 
preserve its topology. As well known from knot theory, local deformations exist on the 
projection of a knot which preserve its topology. These deformations are classified 
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as the Reidemeister moves which are of three types as described in figure ^. 
These topological deformations are sufficient to pass from any two configurations of 
the same knot type. The set of allowed moves in our lattice model will be described 
in detail in the next section. For all the accessible conformations of PSAP of a given 
knot type, we moreover weigh the configuration by a factor 

. . n Q K n w c , x 

7T(P) 



G T (Q,K,w) 

where n is the number of bonds of the PSAP configuration V, with a fugacity K 
per bond, and c is the number of crossings, with a fugacity w per crossing. We also 
artificially introduce a factor rfi with a factor Q > 1 for numerical convenience. This 
factor is unimportant for averages at fixed n (canonical ensemble) but it can improve 
the statistics in a grand-canonical ensemble with varying n. We typically take Q = 2 
in the following. The denominator in eq.(^) is the grand-canonical partition function 
at fixed knot-topology r, which normalizes the weights so that the total weight for all 
accessible configurations is equal to unity. For the knot-topology preserving moves, 
the above weight will dictate the probability to accept or reject the deformation. 



2 Local Elementary Moves 

In this section, we describe the local part of our algorithm, which is a grand- canonical 
implementation, since it involves changes in the number of bonds of the projected 
polygon. Our algorithm uses local elementary moves of four different types depending 
on how many elementary squares of Z 2 (plaquettes) are involved in the move. We 
will use single, double, 3-plaquette (or corner) and 6-plaquette moves. We first give 
here a description of all these elementary moves. The way we decide in the algorithm 
which type of move we attempt and with which probability we accept or reject the 
attempted move is described in the next section. 

Single Moves Single Moves are performed on a single plaquette, and may involve 
length changes through the addition or subtraction of bonds in the PSAP. This 
moves are defined as follows: given a horizontal or vertical bond (pt,pi + i) of the 
current PSAP, we pick a unit vector ti perpendicular to (pi,p i+ i). This defines 
a plaquette (pi,p i+ i,p i+ i + ti,pi + ti). A move can occur only if the edge {pi + 
ti, p i+ i +ej) is not occupied. The move depends on how many edges are occupied 
on the plaquette. If exactly two edges are occupied, we exchange occupied 
and unoccupied edges. This corresponds to a on bead flip (transformation 
Sj in figure ||). If only one edge is occupied, we shift ipi,Pi+\) one lattice 
spacing in the direction of and complete the PSAP by two additional bonds, 
either on the sides, creating a kink (transformation Su in figure |3|), or on 
the diagonals, creating a crossing (transformation Siv or S' IV in figure |J). If 
exactly three edges are occupied, we move (pi,Pi+i) to (pi + ei,p i+ i + ej) and 
delete the two other bonds. The result is either a kink deletion (transformation 
Sin in figure |3[) or a crossing deletion (transformation Sy or S v in figure 
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Figure 3: Single local c-preserving moves. (Si) 
insertion (An = +2). (Sin) Kink deletion (An 



One-bead flip (An = 0). (S n ) Kink 
= -2). 




Figure 4: Single local c-changing moves: Reidemeister I moves. (Siv,S' IV ) Crossing 
insertion (Ac = +1). (Sv,S' v ) Crossing deletion (Ac = —1). 



The single moves can be classified into two different sub-groups depending on 
whether they do or do not preserve the number of crossings c of the PSAP. 
The c-preserving moves illustrated in figure |3] are known in the literature as the 
BFACF moves [|l~3l |i~4| , |j~5[ . While the move Si (one-bead flip) is n-preserving, 
the moves Sn (kink insertion) and Sjjj (kink deletion) change the length of 
the PSAP respectively by An = +2 and An = —2. The c-changing moves are 
illustrated in figure In this case, in addition to the change An = ±2 in the 
length of the PSAP, the number of crossings changes by an amount Ac = +1 
(moves Sjy and S' IV ) or Ac = —1 (moves Sy and S' v ). The moves S' IV and S' v 
are identical to the moves Siy and Sy apart from the reversing of the crossing 
between iS 4 and 1S4. All these c-changing moves correspond to a Reidemeister I 
transformation. 

Double Moves A double move is performed on two adjacent plaquettes. Double 
moves are selected by first choosing, along the current PSAP, a vertex p^. Dif- 
ferent kinds of double moves are considered, depending on the relative orienta- 
tions of the two bonds (pj_i,pj) and (pi,p i+ i) shared by the vertex pi (see figure 

1): 
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Figure 5: Different cases considered for double moves: (a) parallel vertical or hori- 
zontal bonds; (b) perpendicular diagonal bonds; (c) bonds at 135°. 



Case of parallel vertical or horizontal bonds: If the two bonds are par- 



allel, we choose one of the two possible unit vectors perpendicular to 
(Pi-i,Pi+i), and we check if the vertex pj = p^ + belongs to the current 
PSAP configuration. If it does, we look for a 2-plaquette configuration 
T>i such as that depicted in figure If this 2-plaquette configuration is 
encountered, we make the local transformation Dj by exchanging pi and 
Pj = Pi + e ii leading to the configuration T> 2 with two more crossings, or 
make the similar transformation D\ leading to the reversed configuration 



Case of perpendicular diagonal bonds: If {pi-i,Pi) and (pi,Pi+i) are on 
perpendicular diagonal edges, the PSAP necessarily has two crossings in- 
volving two strands of polygon. We check that we have a configuration of 
type T>2 or V 2 i.e. that the two strands are not entangled (see figure^). If 
so, we then make the local transformation Dn or D' u to suppress the two 
crossings. 

Case of bonds at 135°: If one of the bonds, say {pi-i,Pi) is vertical or hor- 
izontal, and the other bond (pi,Pi+i) is at 135° on a diagonal edge, the 
procedure is analogous to the case of parallel bonds, with the unit vec- 
tor perpendicular to (pi-i,pi) and inside the convex sector (pi^i,pi,p i+ i). 
After checking that the vertex pj = Pi + Ci belongs to the current PSAP, 
we look for a 2-plaquette configuration such as the one depicted in figure [7[ 
We then let the crossing diffuse one step, according to transformation Dm 
of figure |7|. A similar procedure is used for the configuration V 3 obtained 
by reversing the crossing of T> 3 . 

Finally, no double move is performed if the two consecutive bonds shared by p^ 
belong to the same plaquette. In this case a corner move (see below) will be 
attempted instead. 

Transformations Dj and Djj can be though of as Reidemeister II moves while 
transformation Dm is simply a diffusion of the crossing to a neighboring pla- 
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Figure 6: Double c-changing moves : Reidemeister II moves. 
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Figure 7: Double c-preserving move: diffusion of a crossing. 



quette. 

In addition to the configurations described above, there are also 2-plaquette 
configurations in which the two involved strands of the PSAP are consecutive. 
We call these configurations degenerate since can be transformed one into an- 
other by the same set of double moves above. Some of such configurations are 
illustrated in figure ||. 

Corner Moves . In these moves, three plaquettes of the underlying squared lattice 
are involved. As for the double moves, we first choose a vertex p^ of the current 
PSAP. A corner move can be performed if the bonds (j)i-x,Pi) and (pi,Pi+i), 
shared by the chosen vertex p i: are on the two consecutive edges of a square, 
as illustrated in figure [| In this case, we check if the vertex pj = —pi + 
Pi+i +Pi-i- belongs to the current PSAP. If so, we look for local configurations 
involving three plaquettes such as those depicted in figure [TI] (Ci,C 2 ) and in 
figure [11] (C 3 ). If the local 3-plaquette configuration does not correspond to any 
of such configurations no move is performed. Otherwise, we attempt one of 
the corner moves illustrated in figures 10 and O. We can think of the corner 
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Figure 8: Degenerate configurations for double moves. 
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Figure 9: Starting configuration for a corner move. 
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Figure 10: Corner c-changing move: Reidemeister II move around a corner. 



moves Cj and Cn as Reidemeister II moves, whereas the corner move of figure 
|TT1 corresponds to a diffusion of a crossing through the corner. Note that this 
diffusion around a corner also involves a rotation of the the crossing on the 
square lattice. Corner moves C\ and C' n connecting configuration C\ to C 2 are 
also possible. In addition, as for the double moves, there are again degenerate 
configurations that can be transformed one into another by the same corner 
moves described above. Some of them are illustrated in figure |T2[ 

6-plaquette Moves (Reidemeister III moves) . This moves are performed by choos- 
ing a diagonal bond in the current PSAP. This bond singles out a crossing. 
We then look for the presence of two neighboring crossings along two adjacent 
diagonal directions, chosen at random. We finally look for configurations such 



as the one depicted in figure [L3[ The performed move Rj corresponds to a 
Reidemeister III move. Here again degenerate configurations exist which can 
be changed one into another. 

Note that double, corner and 6-plaquette moves do not change the number n of bonds 
in the PSAP. 



Pi Pi 




C 3 

Figure 11: Corner c-preserving move : diffusion of a crossing around a corner. 
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Figure 12: Degenerate cases of corner moves. 
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Figure 13: 6-plaquette c-preserving move : Reidemeister III move. 



To end this section let us mention that: 

• For all the moves drawn above, equivalent moves exist which are obtained by 
90° rotations or mirror reflections. 

• If after our checks we do not find one of the above described environments, no 
move is performed. 

• As described in the next section we have an exhaustive procedure to decide 
which type of move will be attempted. 

All the moves above clearly preserve the knot-topology. Therefore any possible 
deformation will maintain the PSAP inside the set of conformations having the same 
knot type as the conformation one started with. In the following, we will assume 
that the algorithm is ergodic, although a full proof of ergodicity goes far beyond our 
goal. By ergodicity, we mean that starting from a given conformation, any other 
conformation with the same knot-topology can be obtained by a series of successive 
elementary moves. In the absence of crossings, our algorithm reduces to the BFACF 
algorithm jT^, [H], which is known to be ergodic. More precisely, the BFACF 
moves (moves Si, Sn and Sm) allow to deform any configuration of a set of self- 
avoiding paths onto any other configuration with the same topology and, if some 
of the paths are open, with the same positions of the end-points. Our algorithm 
completes the BFACF algorithm with knot-topology preserving moves allowing the 
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three types of Reidemeister moves, as well as the migration of a crossing along the 
PSAP. 

Although we have no rigorous proof of ergodicity for our algorithm, let us present 
here what could be the sketch of a proof. Before we proceed, let us note that we 
can limit ourselves without loss of generality to conformations on which none of 
the transformations Sy, Cu or Djj (which decrease the number of crossings) can 
be performed. If such a transformation can be performed, then we do perform it 
(possibly recursively) so as to eliminate the spurious crossings. We are thus left with 
configurations where the suppression of a crossing requires more than one elementary 
move. We then divide the set of remaining conformations with a fixed knot type 
into subsets of conformations which have precisely the same (topologically speaking) 
knot projection. The topology of such a subset is characterized by a fixed number of 
crossings, and a fixed prescription for the connectivity of these crossings. A proof of 
ergodicity can then be established in two steps: 

• (1) Proving the ergodicity within a subset, i.e. for a fixed topology of the pro- 
jection (i.e. a fixed number of crossings and a fixed set of connections between 
them). 

• (2) Showing the possibility to pass from a particular element (representative) 
of one subset to a particular element of another subset if the two different 
projections represent the same knot. 

Since two projections representing the same knot can be related by a series of Rei- 
demeister moves, it is enough to prove point (2) that representatives of two different 
subsets differing by a Reidemeister move can be related one to the other. A good 
choice of representatives is what could be called expanded configurations, i.e. config- 
urations in which all the crossings are well separated. In a grand-canonical ensemble, 
such configurations always exist, whatever big is the number of crossings in the pro- 
jection. A Reidemeister move will involve at most three crossings in the projection. 



For well separated crossings, the figures P, [15| and [16] show how to perform a Rei 



demeister transformation of type I, II or III respectively by use of our elementary 
moves. For a Reidemeister I move, we can use the BFACF moves to shrink a path 
with neighboring end-points onto the segment joining these end-points, then use the 
move Sjy or Sy to create or destroy a crossing, and use BFACF moves again to 
re-inflate the path. For Reidemeister moves II and III, the idea is to use the BFACF 
moves to deform the paths (with fixed end-points) connecting crossings so as to cre- 
ate channels along which the crossings can migrate. The migration of crossings along 
a channel is performed by use of the elementary moves Dm or Cm (at a corner) 
until the crossings come in contact and an elementary move Dj or Du (Reidemeister 
II), or Rj (Reidemeister III) can be performed. The target expanded configuration 
can then be reached again by creating channels, moving the crossings to their final 
positions, and eventually re-deforming the connecting paths. Note that for all these 
deformations, the moves Cj and Cu have not been used. This is because these moves 
are not independent and can be performed by a sequence of D Iy D H) D in , Cm and 
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Figure 14: The passage between two expanded configurations (a) and (b) related by 
a Reidemeister I transformation. Thin lines indicate an arbitrary (supposedly drawn 
a the square lattice) path joining the crossings. For each step of the transformation, 
we indicate the type of elementary move involved. 





Figure 15: The passage between two expanded configurations (a) and (b) related by 
a Reidemeister II transformation. 



BFACF moves. Still in a Monte Carlo program, it is worth implementing them to 
accelerate the algorithm. 

Let us now come to the point (1) of our argument, i.e the ergodicity of the algo- 
rithm inside the subset of conformations having the same projection. A topological 
deformation will involve now only a displacement of the moves and a deformation 
of the paths connecting them, keeping the global connectivity structure. No Reide- 
meister move and in particular no creation of crossings have to be performed. Since 
BFACF moves allow to freely deform the already existing (i.e. with non-zero ini- 
tial length) paths, and the moves Dm and Cm allow to move any crossing along a 
channel of connecting paths, the only possible source of blocking is due to situations 
where crossings have zero length connecting paths, i.e. are nearest neighbors of other 
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Figure 16: The passage between two expanded configurations (a) and (b) related by 
a Reidemeister III transformation. 




Figure 17: A blocked configuration with a loop of 14 alternating crossings. 



crossings, without path between them. Indeed, our moves do not allow us to move 
several crossings at once. For this situation to lead to a real blocking, one moreover 
needs a set of neighboring crossings whose centers form a loop on the dual square 
lattice (the centers of the squares), so that each crossing in the loop in blocked on 
both sides (see an example in figure [IT]). Since spurious crossings have already been 
removed, these loop conformations require moreover many alternating crossings (at 
least eight) and correspond to entangled situations with a complex knot topology. 
We do not expect such conformations to belong to the simple knot topologies that 
we study here where we consider knots made only of trefoils and figure eight primary 
knots. 

3 Updating and Reversibility 

We will now discuss the updating of our system, i.e. how we decide which type of 
move to attempt and with which probability we accept the move. These choices are 
dictated by a condition of reversibility of the algorithm and by the nature of the 
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invariant limiting distribution which we want to reach within an ergodicity class, i.e. 
the weie ht tt(P) of eq.flD which we assign to a configuration V of the PSAP. 

The implementation of the algorithm is "Metropolis style" . Let V\ be the current 
PSAP with ni bonds and c\ crossings. A tentative move of type M is chosen according 
to a procedure which we describe below. It may or may not be possible. If it is, it 
produces a new configuration V2 with length and number of crossings c 2 . Then 
we accept the move with the acceptance rate 



Acc(Ti -> V 2 ) = A(M) min 



hK r 



ri2 

m 



mm 



1, w C2 ~ Cl 



(2) 



The two min[-] terms correspond to the usual Metropolis acceptance rate dictated 
by the limiting distribution of eq. (|]). In particular, the term involving the ratio 
n-ijux is necessary to reproduce the rfi dependence of 7r(P). In addition, we also 
introduce a coefficient A(M) < 1 depending on the type of move M attempted, and 
which we will determine later. The probability to pass from the configuration V\ to 
the configuration P 2 is thus 



Prob(Pi -> V 2 ) = A{M) min 



1 K n2 ~ ni [ — 



mm 



l } w C2 ~ Cl 



xB{V u M) (3) 



where B(Vi, M) is the probability to have selected the move M among all the tentative 
moves which could have been tried on V\. To be sure that the Markov chain converges 
to the limiting distribution 7r (V), we impose on Prob(Pi — ► V2) the detailed balance 
condition 

7T (Vi) Prob{P! -> V 2 } = 7T (V 2 ) Prob{P 2 -> Px} (4) 
This will dictate the choice of A(M) below. 



The implementation of the algorithm is as follows: 



1. To attempt a local move on the current configuration V\ we first select with 
probability q a bond at random on the configuration or with probability (1 — q) 
a vertex at random on the configuration. The number q is a free parameter 
which we can adjust in our simulation. 

2. Case of bond selection. In this case, after having chosen the bond at random 
among the n\ bonds of V\, we attempt a 6-plaquettes move if the bond is a 
diagonal bond, or a single move if the bond is vertical or horizontal. 

(i) 6-plaquette move. Starting from the chosen diagonal bond, we look for a 
configuration such as those in figure [13| (or their degenerate companions). 
If such a configuration is not found the move is rejected. Otherwise, the 
outcome is a PSAP V 2 with n 2 = ri\ and c 2 = Ci, and 

Prob{Pi->P2} = A(i*/) B(Vi,Ri) ; B(Vi,Ri) = ^-x^, (5) 
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The factor 2/nj is the probability to have chosen among n\ bonds one 
of the two bonds forming the crossing which is the starting point of the 
Reidemeister III move, and the factor 1/4 is the probability to have selected 
the correct diagonal directions of the two other crossings involved. The 
important point is that every 6-plaquette move has a unique inverse (see 
figure |13|) by design which is also a 6-plaquette move of the same type. 
Detailed balance holds by choosing 



A{Rj 



1 



(6) 



since, by exchanging V\ and V2, and using 712 — ni, we have ProbjT^ 
Vi} = Prob{P! V 2 } and tt(Pi) = n(V 2 ). 

I Single move. For the single moves, we proceed as described in the previous 
section. If the considered plaquette has two occupied edges, we do a bead 
flip Sj. Since the bead flip Si keeps the number of bonds and crossings 
unchanged, the detailed balance follows by noticing that any such move 
has its own inverse by design. We can thus take 



A(Si) = 1. 



(7) 



If the considered plaquette has a single occupied edge, we must choose be- 
tween the three moves corresponding to an addition of two bonds (Sn,Siv 
and S' IV ) with n 2 — n\ = 2 and Ac = or 1. We choose at random one 
of these three moves with even probability 1/3. For convenience, the ac- 
ceptance rates A(Sn), A(Siv) and A(S' IV ) will be taken all equal to A +2 . 
To fix them, we need to consider the reverse moves, which are performed 
on a plaquette which has three occupied edges. In this case, the current 
configuration dictates without ambiguity which transformation Sm, Sy 
or S' v can be done. Again we take A(Sm), A(Sy) and A(S' V ) all equal to 
A^. If n andP 2 are the two configurations mutually exchanged by some 
of these transformations, with n 2 = n\ + 2, we have 
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the detailed balance 

^i^ViProb-fPi -»■ V 2 } = n$ K n2 w C2 PTob{V 2 -> V x } (11) 
is then satisfied. 
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3. Case of vertex selection. We now suppose that we have selected a vertex, which 
occurs with probability (1 — q). In this case, we attempt either a double or a 
corner move, depending on the relative position of the two bonds (pj_i,pj) and 
iPiiPi+x) shared by the chosen vertex. For the relative orientations of figure [j[ a 
double move is attempted; for the relative orientation of figure ^, a corner move 
is attempted. In both cases the number of bonds remains constant (An = 0). 
For the transformations Dm and Cm in which the number of crossings also 
remains unchanged (figure [7| and figure [11]) the detailed balance holds clearly 
with 

A(Dhi) = A(Cj n ) = 1 (12) 

by noticing that any such move has its own inverse by design. For the remaining 
cases (figure |6| and figure [10]) the number of crossings between V\ = T>i (or 
respectively Ci) and V 2 = £*2 or V 2 (or respectively C 2 or C 2 ) varies by an 
amount c 2 — C\ = 2 and we simply have to account for the fact that there are 
two ways to add the crossings while there is just one way to suppress them. 



Setting A(Z)j) = A(d) = A%' C) and 2A(D H ) = A(C n ) = A^>, we have 



{D,c) ■ 



Prob{P! — > V 2 } = A 
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x - x mm 
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1,W 



(13) 



where the first line is for a double move and the second line for a corner move. 
The extra coefficients 2 and 1/2 in the case of a double move come respec- 
tively from the two different vertices which select the pair of parallel strands 
under consideration and from the probability 1/2 to make the correct choice of 
direction e; in which we move the vertex pj. We have conversely 



Prob{P 2 -> Vi} 
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The detailed balance 



n Q K nl w Cl Prob{Vx -> P 2 } = n$ K n2 w C2 PTob{V 2 
is then satisfied by choosing 
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(15) 



(16) 



The analysis above may seem a little cumbersome. In practice, the only subtle 
point is the determination of the coefficients A(M), the list of which is displayed 
in Table 0. These coefficients simply correct the fact that some transformations 
are in competition (like a kink creation and a crossing creation) while the reverse 
transformation has no competitor. As usual, the Metropolis criterion simply creates 
the correct "energy" factor. With the above choices for the acceptance coefficients 
A(M), the transitions probabilities verify the detailed balance condition (|J) for all 
PSAP's Vi, V 2 , with the invariant limit distribution of eq.(|l]) 
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Single Moves 


Double Moves 


Corner Moves 


6-plaquette Moves 


Si 


1 


Di 


1 


Ci 


1 


Ri 


l 


Sii 


1 


D u 


1/4 


C n 


1/2 






Shi 


1/3 


Dm 


1 


Cm 


1 






Siv, S' IV 


1 














Sy, S'y 


1/3 















Table 1: Coefficients A(M) for the different local moves. 



4 Improvements of the algorithm 
4.1 Non-local Moves 

Algorithms based on local moves are known to have very long autocorrelation times 
which in turn lead to large statistical errors. For instance, it is known that, for the 
BFACF algorithm (which in our language simply corresponds to forbid crossings), 
the so called "exponential" autocorrelation time, which controls the relaxation of the 



Markov chain from an initial configuration to equilibrium, is infinite jig] . For the 
same algorithm, the so called "integrated" autocorrelation time, which controls the 
statistical error in measured quantities, is usually very large. 

To remedy this problem, algorithms have been proposed [Q, [HJ which combine 
local moves and non-local moves. The non-local moves correspond to select two 
points far apart on the chain and to reverse the part of the chain between them. In 
these hybrid algorithms, the non-local moves hopefully ensure the rapid equilibration 
within subspaces of fixed length n, while local moves ensure equilibration between 
different values of n. Now if one wants to study the statistics of polygons with a fixed 
knot type, one must make sure that the non-local moves preserve the topology of the 
polygon. For the usual non-local moves performed on 3-dimensional Self-Avoiding 
Polygons in Z 3 , this is in general not the case. One way to deal with this problem is 
then to check potential changes in the knot type by calculating a topological invariant 
such as the Alexander polynomial at every successful non-local move 0. This is 
however computationally costly and it also has the disadvantage that the Alexander 
polynomial is not a perfect characterization of the knot type and distinct knots can 
have the same Alexander polynomial. 

The situation is different in our model where the two-dimensional character of the 
PSAP allows the introduction of non-local moves that do preserve the topology of the 
PSAP. The non-local move we consider are performed by first selecting two vertices 
(pi and P2) in the PSAP with uniform probability. If these vertices are neighbors, no 
move is attempted; if not, the vertices separate the PSAP into two arcs, say a± and 
ci2- If the two arcs have no crossings in common we can perform a rotation by 180° of 
the shortest one, say a\ (see figure |18|). Let a[ be the image of a± under this rotation. 
If a[ and 02 do not intersect each other, no new crossing has been generated and the 
move is accepted. The important feature of this move is that it does not change the 
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Figure 18: Non-local transformation of a PSAP. The arc a[ is the image of ai by a 
rotation of 180°. The arc a 2 is unchanged. 






Figure 19: The passage from a trefoil to its image by a rotation of 180° by use if local 
Reidemeister moves only. We indicate only the major intermediate steps. 




Figure 20: The passage from a figure eight to its image by a rotation of 180° by use 
if local Reidemeister moves only. We indicate only the major intermediate steps. 
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topology of the arc a\ and therefore, since we avoid topological interaction between 
the arcs, it preserves the topology of the entire PSAP. To be more precise, the absence 
of crossings between the two arcs guarantees that the knots in a x and those in a 2 are 
well separated and the total knot can be viewed as the compositions of the two knots. 
It is known that the composition of two knots can give at most two different results 
according to whether their (arbitrarily fixed) orientations do or do not match in the 
composition. Note that our transformation reverses this matching of orientations. 
The result is however unique as soon as at least one of the knots is reversible (i.e. its 
two opposite orientations are in fact equivalent). This is the case for all the knots 
that we will use, which are made of trefoils (3i) and figure eight (4i) knots only, both 
of which are reversible. In practice, we can also argue that, for a trefoil (3i) or figure 
eight (4 X ) knot, the passage from the arc a x to its image a[ can be done by using 
Reidemeister moves which do not involve the arc a 2 . This is illustrated for the trefoil 
knot in figure [B1 and for the figure eight knot in figure 120. 



The non-local move described above is the simplest one that can be implemented. 
More sophisticated moves can be considered in which the constraint of having no 
crossings between a\ and a 2 can be relaxed to weaker conditions. We however limited 
ourself to the simplest non-local moves without crossings. 

In our algorithm, we therefore have an additional free parameter < p < 1 which 
is the probability that we attempt non-local move rather than a local move. 

4.2 Multiple Markov Chain Method 

Another way to improve the slow convergence of the local algorithm is by sampling 
along a set of several Markov chains run in parallel. 

This method, refered to as the Multiple Markov Chain (MMC) method, has been 
introduced by Geyer JTB| and then adapted to interacting self-avoiding walks with 



fixed length by Tesi et al |19[ (see also J20J for a more detailed analysis of the method). 

We describe here a variant of this approach suitable for simulations in the grand 
canonical (n varying) ensemble. The starting point is the single Markov chain whose 
limiting distribution 7r(V) depends on the step and crossing fugacities K and w. 
Here we will consider a procedure with a fixed value of w, hence the w dependence 
will be omitted in this section. For values of K close to (and below) the critical 
value K c where the average length of the PSAP diverges, the Markov chain samples 
configurations with very large values of n and the convergence of the algorithm is 
therefore very slow. On the contrary for a value K' of the step fugacity such that 
K' « K c , the sampled configurations have small values of n and the convergence 
of the local algorithm is fast. The idea is to select a set of values K' = K\ < K 2 < 
K 3 < . . . < K m = K to interpolate between K' and K, so that Kj and Kj + i are close 
enough to ensure that there is a considerable overlap between the distribution at Kj 
and that at Kj+i. 

The m Markov chains at K±, K 2 , . . . , K m are evolved in parallel, and we let the 
chains interact by possibly exchanging conformations as follows: select at random 
two neighboring chains with fugacities Kj and Kj + \. A trial move is an attempt 
to swap the two current conformations of these chains. If we denote by itj(V) the 
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K 


(1) local 


(2) local + non-local 


(3) local + MMC 


(4) all 


(n) 


r(n) 


(n) 


r(n) 


(n) 


r(n) 


(n) 


r(n) 


0.3760 


76 ±4 


18 ±3 


75 ±3 


12 ±2 


75 ±1 


2.2 ±0.2 


75 ±1 


1.5 ±0.2 


0.3771 


110 ±11 


50 ± 10 


100 ±6 


20 ±5 


111 ±4 


5.6 ±0.9 


112 ±2 


3.0 ±0.3 


0.3780 


160 ± 19 


67 ±11 


188 ± 18 


40 ± 15 


188 ±9 


14 ±3 


192 ±9 


13 ±3 


0.3783 






260 ± 35 


90 ±25 


231 ± 17 


47 ± 15 


259 ± 19 


35 ± 10 



Table 2: Integrated autocorrelation times for the average length (n) of the unknot 
at w — 1, for (1) the purely local algorithm, (2) the local + non-local algorithm, (3) 
the local + MMC algorithm and (4) local + non-local + MMC algorithm. Each run 
has a total of 10 8 local moves and the unit of sampling is 10 4 local moves. The data 
which are not reliable are omitted. 



equilibrium weight of the state V for the chain at fugacity Kj, and Vj and Vj+i the 
current conformations in the j-th and (j + l)-th chain, we accept the trial move (i.e. 
swap Vj and Vj+i) with the acceptance rate 

a (V^V 1+1 )=^L^^^\. (17) 

Note that we do not need to compute the normalization of the limiting distribu- 
tions (i.e. the partition functions) since these normalizations drop out in the ratio 
above. 

The whole process is itself a Markov chain, which we can call a composite Markov 
chain. Since the underlying Markov chains are ergodic (for a fixed knot type), so is the 
composite Markov chain. Moreover the composite chain is in detailed balance since 
the "swap" -move as well as the moves in the underlying chains are. Consequently, 
the invariant limiting distribution is the product of the distributions of m separate 
Markov chains at the temperatures K i < K 2 < . . . < K m . 

One can clearly understand the advantage of this method by focusing on a par- 
ticular value Kj. Every successful swap involving this chain and the chain Kj+i 
corresponds to a big change in the configuration at this value Kj, so the correlation 
time at each Kj is markedly reduced. Swaps move the system at a particular Kj to 
new regions of the configuration space and therefore tend to eliminate problems of 
quasi-ergodicity. If two chains Kj and Kj + \ are far apart, the changes in the config- 
uration are likely to be large but eq.fll7|) suggests that bigger distances between the 



two selected chains lead to smaller acceptance rates for the swaps. In practice, the 
set of the values of K is chosen so that successive values of K are sufficiently close for 
swapping to occur rather frequently, and K\ is small enough for convergence at K\ 
to be rapid. All this requires that the total number m of chains is sufficiently large. 
Hence the gain in convergence rate is obtained at only marginal cost in computer 
time. 

In Table |2| we compare our estimates of the integrated autocorrelation time (given 
in units of sampling, i.e. per 10 4 attempted local moves in our tests) of the average 
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PSAP lengths, as a function of K, for the four different implementations of the 
algorithm: (1) the purely local algorithm using only elementary local moves, (2) the 
local + non-local algorithm, using both local and non-local moves, (3) the local MMC 
algorithm, using local moves only together with the MMC sampling and (4) the local 
+ non local MMC algorithm, using local and non local moves together with the MMC 
sampling. In all cases, the integrated autocorrelation time for the polygon average 
length increases as K approaches the critical value K c , as expected. However, it is 
clear that both non-local moves and MMC swaps lead to a significant decrease in this 
integrated autocorrelation time. 

5 Results for w = 1 

In this section, we focus on the model with w = 1, i.e. we do not attach any particular 
weight to the crossings. A typical configuration with the topology of a trefoil knot is 
shown in figure [H]. 




Figure 21: A typical configuration for the 3i at w — 1. The number of bonds is 
n = 906. We have indicated the location of the knotted section. 



5.1 Connectivity /i and Entropic Exponent a 

We present here our estimates for the connectivity constant \x of the PSAP and the 
entropic exponent a for w = 1. These quantities are defined by the large n behavior 
of the number p n {r) of configurations of fixed knot-topology r with a fixed number 
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n of bonds, expected to be of the form 

Pn(r) ~ A(t) (/i(r))" n"M~ 3 (18) 

In particular we study the dependence of /x(r) and a(r) on the topology r of the 
PSAP. A method to compute /x(r) and «(t) relies on the asymptotic behavior for the 
the average length (n(r)) of the PSAP for K close to the critical fugacity K c (t) = 
(/x(r)) -1 , namely: 

This behavior directly follows our choice of the statistical weight in eq.(|T]), which 
yields for w — 1 

<n(T)> = E„^.,(r)A- (20) 

Inserting eq.([TBJ) directly leads to eq. (IDj). At leading order, we can use eq. (|T9"D to 
approximate l/(n(r)) as: 

1 ~ 1-^) _ 1 1 (21) 



(n(r)> [a(r) + Q-2]^(r) [a(r) + Q- 2]Kh(t) a(f) + Q-2' 

An estimate of /i(r) and a(r) can therefore be obtained by a linear extrapolation of 
l/(n(r)) as a function of for l/(n(r)) — > 0. 

In figure we have plotted 1/ (n(r)) as a function oil/K for the topology of the 
unknot (0), that of the trefoil knot (3i), that of the figure eight knot (4i), that of the 
composite knot made of two trefoils (3i#3i) and that of the composite knot made of 
a trefoil and a figure eight (3i#4i). These data are obtained with Q = 2. We note 
that as the complexity of the topology increases, the corrections to the linear scaling 



(p|) become more and more important, but for K sufficiently close to the critical 
value K c (t) = /i(r) _1 , a linear behavior is indeed obtained. By extrapolating the 
data, using the linear fit of equation fl21~|), we obtain the following estimates: 



/x(0) = 3.254 ±0.01, 

/x(3i) = 3.250 ±0.06, 

/i(4i) = 3.255 ±0.07, (22) 

fz(3 x #3i) = 3.249 ±0.13, 

/x(3i#4 x ) = 3.261 ±0.15. 

These values coincide up to the second decimal, and it seems reasonable to assume 
that they are indeed all equal. In fact, since they were measured from completely 
independent simulations, we can take their average to estimate the growth constant 
of PSAP of a fixed knot type: we obtain /x(r) = fii = 3.251 ± 0.002 (95% confidence 
interval) for PSAP's with arbitrary fixed topology r. From the same linear fit an 
estimate of the entropic exponent a can also be obtained: 

a(0) = 0.58 ±0.07, 
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2.10 ±0.08, 

2.13 ±0.10, (23) 
3.91 ±0.24, 
4.05 ±0.25. 

In contrast with /i, the exponent a depends clearly on the topology. However, this 
dependence appears to be weak in the sense that a depends only on how many primary 
knots form the knot under consideration and not on the precise nature of these knots, 
hence a(3i) = a(4 x ) and a(3i#3i) = a(3i#4i) within numerical uncertainties. Such 
a result was already found in [|1(], [Tl|] for three-dimensional knotted polygons. There 
however, it was found that a(r) increases by one for each added primary knot, leading 
to ol(t) = ot(0) + Nf where r is a knot with Nf prime factors. This formula does not 
seem to hold in our case. 



a(3i) = 

a(4i) = 

a(3i#30 = 

a(3i#4i) = 
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Figure 22: Plot of 1/ (n) as a function of l/K for the unknot (•), the trefoil (o), the 
figure eight (v), the composite knot 3i#3i (a) and the composite knot 3i#4i (#). 



Another way to estimate fi and a relies on the maximum likelihood method as 



suggested in [pi]] . In this case, we assume that for n > n min the relation 

p n = ^ n «- 3 (1 + a /n) (24) 

holds, where a is a free parameter that we can vary to mimic a correction to scaling. 
By varying n m i n and a we can probe the sensitivity of our estimates to scaling. We 
generally find a region of minimum sensitivity from which we extract the following 
estimates: 

//(0) = 3.254 ±0.005, 
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/i(3i) = 3.251 ±0.07, 

^(4i) = 3.253 ±0.09, (25) 

//(3i#3i) = 3.245 ±0.10, 

/x(3i#4i) = 3.255 ±0.11. 



and for the entropic exponents: 



a(0) = 0.58 ±0.03, 
a(3i) = 1.86 ±0.06, 



a(4i) = 1.95 ±0.08, (26) 
a(3i#3i) = 4.20 ±0.15, 
a(3i#4i) = 3.90 ±0.20. 

All these results are in agreement with the estimates (E3) and fl2~4]). 



5.2 The Metric Exponent v 

For the mean square radius of gyration, we assume the large n behavior 

(R 2 (r)) n = A v {r)n 2 ^ (27) 

where the average is performed over configurations with fixed n. In order to estimate 
v from our data in the grand-canonical ensemble, we compute the expected value of 
the mean square radius of gyration, 

{R (T)) = LnW (M) 

Using the asymptotic forms ( |i~8"|) and (p7[), we get the following approximation for 
(R 2 (t)) close to the critical K c : 

{R 2 (r)} ~ h(r) (1 - K^(t)Y 2v{t) (29) 

or equivalently, using eq. ([T9|) , the approximation: 

(R 2 (r)} -h'(r)(n(T)} 2u{r) (30) 

where we now have a relation between (R 2 (t)) and (n(r)). A log-log fit of our data to 
equation (|30D , where .fT takes a different value for each chain in our multiple Markov 
chain Monte Carlo, can be used to estimate v for each knot type. Our best estimates 
are 

i/(0) = 0.747 ±0.005, 
u(3x) = 0.758 ±0.010, 

i/(4i) = 0.748 ± 0.010, (31) 
i/(3i#3i) = 0.760 ±0.020, 
i/(3i#4 x ) = 0.720 ±0.040. 
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Within the stated error bars, the results above are all identical and consistent with the 
exact value z^saw = 3/4 |23j for self-avoiding walks or self-avoiding polygons in d = 2. 
This seems to indicate that although polygons in our model are not strictly self- 
avoiding since crossings are allowed for w = 1, the presence of these crossings limited 
to diagonal edges does not modify the scaling behavior of the radius of gyration. 

6 Results for w « 1 

It is interesting to consider our model in the limit w — > where crossings become 
forbidden. At w — 0, our model coincides with that of a strictly self-avoiding polygon 
on the square lattice. Still, we can start from a configuration with a non-trivial knot 
topology r (say that of the trefoil knot) with its minimal number c m i n (r) of crossings 
(c m in(3i) = 3 for the trefoil) and look at the limit w — > of G T (Q, K,w)/w Cmin( - T \ In 
this case, we explore the phase space of all the configuration of knot type r and with 
the minimal number of crossings c m i n (r). This limit is the "minimal" extension of 
the strictly self-avoiding model to include a non-trivial topology. In this limit, we of 
course expect the model to be in the universality class of strictly self-avoiding walks. 

6.1 Results for fi and a 

In practice, the easiest way to study the w —>■ limit is simply to fix in our simulation 
a very small value of w. The results presented here have been obtained with the value 
of the crossing fugacity w fixed to w — 0.0001. Here again we define the connectivity 
constant fi and the entropic exponent a for the asymptotic behavior of the number 
Pn\ T ) °f configurations of polygons in the topology r, with n bonds and with the 
minimal number c m i n (r) of crossings compatible with r. As for the case w = 1, we 
can estimate \x and a by plotting l/(n(r)) as a function of 1/K and by making a 
linear extrapolation for 1/ (n(r)) — > 0. 

For fi we obtain the following estimates: 

2.640 ±0.01, 

2.640 ±0.01, 

2.639 ±0.02, (32) 

2.641 ±0.05, 
2.645 ±0.08. 

These values coincide to the second decimal place, and it seems reasonable to assume 
that they are indeed all equal. It is interesting to notice that these estimates of fj,(r) 
are a good upper bound for the best available estimates of the connectivity constant 
for self-avoiding polygons in the d — 2 square lattice (/x = 2.63815) p2"| . 
For the entropic exponents we obtain: 

a(0) = 0.503 ±0.040, 
a(3i) = 2.24 ±0.10, 



M0) = 

M3i) = 

M4i) = 

M3i#3i) = 

M3i#4i) = 
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a(4i) = 2.28 ±0.10, (33) 
a(3i#3i) = 3.2 ±0.3, 
a(3i#4i) = 3.3 ±0.5. 

As for the case w — 1, we again see a dependence of a on the number of prime 
factors forming the knot. We note however that the values of a for a fixed knot type 
are different in this limit (w — > 0) from the w = 1 values. For the unknot, we recover 
the exact value «saw = 1/2 j2"3"| . This is expected since strictly self-avoiding rings 
in two dimensions cannot form knots, thus the usual a exponent derived form field 
theory corresponds to fixing the topology to that of the unknot rather than summing 
over all topologies as in three dimensions. 



6.2 The Metric Exponent v 

As for the case w — 1, we have evaluated the exponent v in the limit w —>■ from 
the relation of eq. ^between the average radius of gyration (R 2 (r)) and the average 
number of bonds (n(r)}. Our estimates are 

i/(0) = 0.748 ±0.015, 
z/(3i) = 0.753 ±0.010, 

z/(4i) = 0.749 ±0.011, (34) 
i/(3i#3i) = 0.736 ±0.015, 
i/(3i#4 x ) = 0.742 ±0.020. 

We thus obtain an exponent v independent of the knot-type and in perfect agree- 
ment with the exact value t'sAw = 3/4 for self- avoiding walks in two dimensions. 



6.3 Statistics of Arc Lengths: Localization of Knots. 

One advantage of taking this limit w — > is that, by preventing the creation of 
spurious crossings, we have a fixed minimal number c m ; n (r) of crossings which we can 
easily localize along the PSAP. In particular, we can study the typical "size" of the 
knot in a given topology and see for instance if primary knots tend to be localized 
on a small number of bonds or prefer to extend over the whole PSAP. A simple way 
to analyze this phenomenon is to study the distribution on the lengths of the arcs 
between consecutive crossings. 

In the case of a trefoil knot, we have exactly c m j n (3x) = 3 crossings which separate 
exactly 6 arcs along the PSAP. For each configuration with n bonds, we can classify 
these arcs in order of decreasing arc length U = Xin with Ai > A2 > ■ • ■ > A6 and 
Yh=i Aj = 1. Here A, is thus the fraction of the total length n inside the i-th largest 
arc. 

In figure E3L we have plotted the fraction of length occupied by the two largest 



arcs for a trefoil knot, as a function of K. We see that as K — > K c ((n) — > 00), Ai — > 1 
while A2, and thus all the other Aj for % > 2, tend to zero. This result is consistent 
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Figure 23: Average length fraction Ai and A2 of the two largest arcs along the PSAP 
as a function of K for the trefoil (Ai — > o, A2 — > □), the composite knot 3i#3i (v and 
a) and the composite knot 3i#4i (• and *). The dashed lines indicate the values 3/4 
and 1/4 expected for two strictly localized knots moving freely along the PSAP. 



with a localized knot which takes place over a number a bonds which is either finite, 
or at least growth less rapidly than n. This result of localization is corroborated by 
a direct visualization of a typical configuration, as shown in figure p4] . 

Another interesting case is that of composite knots 3i#3i and 3i#4 1; which have 
respectively c min (3i#3i) = 6 and c min (3i#4 1 ) = 7. These knots are thus made 
respectively of exactly 12 and 14 arcs. We have plotted on figure 23 the proportions 
of the total length n occupied by the two largest arcs in both cases. Although the 
statistical errors are rather important when K — > K c , we see in each case that these 
proportions add up to one, with the proportion occupied by all the smaller arcs 
tending to zero. This result is consistent with a picture of the composite knot as 
made of two perfectly localized primary factors. The localization of the two primary 
knots is also directly visible on a typical configuration of the PSAP (see figure |25|) . 
If we assume that these localized factors move freely along the PSAP and interact 
only locally by preventing to pass through one another (which would require creating 
additional crossings), we should expect that the proportions tend to 3/4 and 1/4. 
Indeed, for two points at random drawn on a loop of size n, the largest arc between 
them has average length (1/n) Jq dx max(x, n — x) = (3n/4) and the complementary 
arc has thus average length n/4. In practice, we observe a value a Ai < 3/4 and 
A2 > 1/4, which tends to indicate that the primary knots repel each other. This 
repulsion might however be a finite size effect. Indeed the deviation from the expected 
asymptotic values 3/4 and 1/4 is larger for the 3i#4i knot than for the 3i#3i knot 
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Figure 24: A typical configuration for the trefoil knot (3i) at w ~ 0. The number of 
bonds is n = 312. We have indicated the location of the knotted section. 



as the Ai prime factor typically requires more bonds to be formed. 

7 Results for w > 1 

Up to now, we have studied only the case w — 1 where no fugacity is associated to the 
number of crossings and the case w — > where we have a minimal number of crossings 
compatible with the knot topology at hand. Although we obtained different values 
for the (non universal) connectivity constant //, we found in both cases a value of v 
compatible with the exact result t'sAw = 3/4 in two dimensions, which suggest that 
these two cases are in the universality class of self-avoiding walks, and so is probably 
any case of the model with a value of w between and 1. In this section, we will 
rather be interested in the behavior of the model for w > 1. 

7.1 Connectivity \i 

Let us first concentrate on the connectivity constant \x. If we denote by p n ,c{ T ) the 
number of configurations of PSAP with n bonds, c crossings and the knot topology 
r, we can introduce the canonical partition function 

Pn (w,T)= £ Pn , c (r)w C - C ^ (35) 

For w — 1, we recover the number p n (r) of configurations with knot-topology r and 
with n bonds, irrespectively of their number of crossings. For w = 0, we recover the 
number p$(r) of configurations with topology r, number of bonds n, and with the 
minimal number of crossings c m i n (r). 

We define n(w, r) and a(w, r) by the large n behavior 

p n (w, t) ~ A(w, t)h(w, T ) n n a(w ' r) - 3 (36) 
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Figure 25: A typical configuration for the composite knot 3i#3i at w ~ 0. We have 
indicated the location of the two primary knots. 



In the grand canonical ensemble, we have 

G T (Q,K,w)/w c ^ T) = J2Pn(w,T)K n n Q 



(37) 



hence we can again relate /j,(w,t) to the inverse of the critical value K c (w) of K for 
which (n(w, r)) diverges. 

Before we present our numerical results, let us first establish a lower and an upper 
bound for [x{w). 

Lower bound for w > 1: The canonical partition function p n (w,r) is clearly 
bounded from below by the contribution w c of any acceptable configuration in the 
given knot-topology class. Since w > 1, we can obtain a larger lower bound by se- 
lecting a configuration with the maximum number of crossings. To obtain such a 
configuration, we start from any particular realization of the knot, as compact as 
possible. It will have, say, n bonds and c crossings. On any outermost bond, we 
can start to grow a linear sequence of crossings, such as the one depicted in figure 



26| . Note that we can arbitrarily choose any of the two possible crossings for each 
new crossing added in linear piece. The addition of this linear sequence clearly does 
not modify the topology of the knot. If n is the total number of bonds in the con- 
figuration, we get a total number of crossings equal to c + {n — n )/2 for any of the 
2( n - n o)/2 guch configurations. We thus get 

l(r) x (2w)i (38) 



Pu{w,t) > 2 2 w 
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Figure 26: The growing of a linear sequence of crossings from a compact knot. 



In the large n limit, this directly leads to the lower bound 

fi(w) > V2w (39) 

Upper bound for w > 1: It is clear that any configuration with n bonds can have 
at most n/2 crossings, irrespectively of its knot type. We thus get for w > 1 

p n {w, t) < p n (r)w% (40) 

This directly leads to 

/ji(w,t) < frVw (41) 



From eqs.(p9|) and (|l|), we can deduce that at large w, n(w) behaves like 

fi(w, r) w ~°° X (42) 

with ^2 < A < ~ 3.251. 

To complete our analysis, we can also give lower and upper bounds on the case 
< w < 1. In this case, it is clear that 

p^(T)<p n (w,T)<p n (r) (43) 

and therefore 

2.638 ~ < n(w, t) < fi x ~ 3.251 (44) 

Figure ^ presents our results for fi{w) together with the above lower and upper 
bounds. At large w, we confirm the behavior of eq.(H^) with A ~ 1.46 close to the 
lower bound v2. 

7.2 The Metric Exponent v\ Branched Polymer Behavior 

We have also estimated the metric exponent v for different values of w > 1. Our best 
estimates are 

i/(u; = 2) = 0.668 ±0.008, 
u(w = 3) = 0.665 ±0.007, 
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Figure 27: Estimate of function of the crossing fugacity w; the solid lines 

correspond respectively to the curves \i = 3.2511/7 2 (upper bound) and /i = (2u>)2 
(lower bound). The dashed line corresponds to the fit /i = 1.46(u;2 — 1) + 2.86. 



v{w - 


= 4) = 


0.645 ±0.008, 


u(w = 


10) = 


0.634 ±0.010, 


v{w = 


15) = 


0.630 ±0.011, 


u{w = 


20) = 


0.625 ±0.012, 


u{w = 


30) = 


0.620 ±0.015, 


v{w = 


50) = 


0.624 ±0.020, 



(45) 



These values are plotted in figure £8|, together with the values estimated for w — > and 
w — 1. Within the error bars, we distinguish two values for v. For w < 1, our data are 
consistent with the value i^saw = 3/4 of self-avoiding walks. For w > 1 however, our 
estimates are far below this value and approximately located around the value z/bp ~ 
0.64 obtained for branched polymers in D = 2 0, ^|. At large w, the branched 
polymer picture is corroborated by looking at a typical configuration, as the one 



depicted in figure |29|. Such configurations are naturally selected at large w since long 
one-dimensional coils maximize the number of crossings and are thus energetically 
favored, while branching points, occuring at a finite energy cost, are favored for 
entropic reasons. From our data, a transition from the self-avoiding behavior to a 
branched polymer type statistics seems to occur at w = 1. For w < 1, crossings 
are penalized and play a marginal role without changing the asymptotic behavior of 
the polymer, which can be considered as self-avoiding. On the contrary, for w > 1, 
a finite density of crossings is present in a typical configuration, which degenerates 
into a branched polymer conformation. At the transition point w — 1, we find 
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Figure 28: The exponent v as a function of w. The dashed lines indicate the values 
of v for self-avoiding walks and branched polymers. 



occasionally interesting configurations showing a phase separation between a phase 
rich in crossings and a phase with few crossings (see figure ^). In this case, the 
largest of the largest value v = vsaw is observed. Finally, let us note that a similar 
transition from a self-avoiding behavior to a branched polymer behavior has already 
been observed in an interacting two-tolerant trail model on the square lattice ||26|| . 

8 Conclusions 

In this paper, we introduced a two-dimensional lattice model of Projected Self- 
Avoiding Polygons, describing two-dimensional knotted polymers rings with a fixed 
knot topology. We made a numerical study of the model using a Monte Carlo algo- 
rithm based on both local and non-local moves which preserve the knot topology of 
the polygon, together with a Multiple Markov Chain procedure. In the absence of a 
rigorous proof of the ergodicity of our algorithm within a fixed knot type, we gave 
strong arguments in favor of this ergodicity, together with the guidelines of a possible 
proof. We presented our numerical results for the connectivity constant /x and the 
critical exponents v and a for various topologies and various values of the fugacity w 
associated to the number of crossings in the projected polygon. The connectivity \i is 
found to depend only on this fugacity w and not on the particular knot topology at 
hand. This is somewhat expected since \x typically measures a local property which 
is not affected by a global constraint such as fixing the knot topology By measuring 
the radius of gyration of the polygon, we observe only two possible values for the 
exponent v: a self- avoiding value 3/4 found at low w < 1 irrespectively of the knot 
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Figure 29: A typical configuration for the unknot (0) at w — 30. 




Figure 30: A configuration for the unknot at w — 1, with n = f250 bonds. One 
clearly sees two domains with very different densities of crossings, one being sparse 
while the other is dense. 
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type, and a large branched polymer value ~ 0.64 found at large w > 1 for all knot 
topologies. Our data suggest that the transition between the two regimes occurs 
precisely at w — 1. Finally, at fixed w, the exponent a for the number of accessible 
configurations is found to depend only on the number of prime factors forming the 
knot. A weak w dependence is also observed. 

One nice feature of our two-dimensional model is the possibility to have more 
specific characterizations of the polygon conformations, such as the number of cross- 
ings or the arc lengths between successive crossings along the chain. This allowed us 
in particular to measure the typical extension of a knot inside the polygon. In the 
limit w — > 0, we found for instance that the prime factors in a knotted polygon are 
localized and weakly interacting. 

Other applications of this model can be imagined, such as a statistics of the knot 
type of randomly drawn polygons. For an arbitrary polygon, the determination of 
the knot topology could be simply achieved in principle by imposing a very low value 
of w and letting the system evolve toward a configuration with a minimal number of 
crossings in the projection. In particular, it would be interesting to understand how 
the knot complexity increases with the number of monomers in the chain. 



34 



References 



[I] P. G. de Gennes, Exponents for the excluded volume problem as derived by the 
Wilson method Phys. Lett. A 38 (1972) 339-340. 

[2] D.W. Sumners and S.G. Whittington, Knots in self-avoiding walks, J. Phys. A: 
Math. Gen. 21 (1988), 1689-1694. 

[3] N. Pippenger, Knots in random walks, Discrete Applied Math. 25 (1989), 273-278. 

[4] V.V. Ribenkov, N.R. Cozzarelli and A.V. Vologodskii, Probability of DNA knotting 
and the effective diameter of the DNA double helix, Proc. Natl. Acad. Sci. U.S.A. 
90 (1993), 5307-5311. 

[5] S.Y. Shaw and J.C. Wang, Knotting of a DNA chain during ring closure, Science 
260 (1993) 533-536. 

[6] E. J. Janse van Rensburg and S.G. Whittington, The dimensions of knotted poly- 
gons, J. Phys. A: Math. Gen. 24 (1991), 3935-3948. 

[7] T. Deguchi and K. Tsurusaki, Topology of closed random polygons, J. Phys. Soc. 
Japan 62 (1993), 1411-1414. 

[8] T. Deguchi and K. Tsurusaki, Statistical study of random knotting using Vasiliev 
invariants, in Random Knotting and Linking ed K C Millett and D. W. Sumners 
(Singapore: World Scientific) (1994), 89-121. 

[9] S. Quake, Topological effects of knots in polymers, Phys. Rev. Lett. 73 (1994), 
3317-3320. 

[10] E. Orlandini, M.C. Tesi, E.J. Janse van Rensburg and S.G. Whittington, En- 
tropic exponents of lattice polygons with specified knot type, 

[II] E. Orlandini, M.C. Tesi, E.J. Janse van Rensburg and S.G. Whittington, Asymp- 
totics of knotted lattice polygons, J. Phys A: Math. Gen. 31 (1998), 5953-5967. 

[12] For an introduction to knot theory, see for instance C.C. Adams, The Knot Book: 
an elementary introduction to the mathematical theory of knots, Ed. W. H. Freeman 
and company, New York (1994). 

[13] B. Berg and D. Foester, Random paths and random surfaces on a digital com- 
puter, Phys. Lett. 106B (1981), 323-326. 

[14] C. Aragao de Carvalho and S. Caracciolo, A new Monte Carlo approach to the 
critical properties of self-avoiding random walks, Journal de Physique 44 (1983), 
323-331. 

[15] C. Aragao de Carvalho, S. Caracciolo and J. Frohlich, Polymers and g\<fi\ A -theory 
in four dimensions, Nucl. Phys. B [FS7] 215 (1983), 209-248. 



35 



[16] A.D. Sokal and L.E. Thomas, Exponential convergence to equilibrium for a class 
of random walk models, J. Stat. Phys. 54 (1989), 797-828. 

[17] S. Caracciolo, A. Pellissetto and A.D. Sokal, Nonlocal Monte Carlo algorithm for 
self-avoiding walks with fixed endpoints, J. Stat. Phys (1990), 1-53. 

[18] C.J. Geyer, Markov chain Monte Carlo maximum likelihood, in Computing Sci- 
ence and Statistics: Proceedings of the 23rd Symposium on the Interface, ed E M 
Keramidas (Interface Foundation) (1991), 156-163. 

[19] M.C. Tesi, E.J. Janse van Rensburg, E. Orlandini and S.G. Whittington, Monte 
Carlo study of the interacting self-avoiding walk model, J. Stat. Phys. 82 (1996), 
155-181. 

[20] E. Orlandini, Monte Carlo study of polymer systems by Multiple Markov Chain 
method, IMA Proc. Work 7 (1996). 

[21] A. Berretti and A. Sokal, , J. Stat. Phys. 40 (1985), 1689-1694. 

[22] S. Caracciolo, A. Pellissetto and A.D. Sokal, Monte Carlo test of hyperscaling 
relation for the two-dimensional self-avoiding walk:II, J. Phys. A: Math. Gen. 23 
(1990), 4509-4517. 

[23] B. Nienhuis, Exact critical point and critical exponents of O(n) models in two 
dimensions Phys. Rev. Lett. 49 (1982), 1062-1065. 

[24] B. Derrida and L. De Seze, Application of the phenomenological renormalization 
to percolation and lattice animals in dimension 2, J. Physique 43 (1982), 475-483. 

[25] E.J. Janse van Rensburg and N. Madras, A non-local Monte Carlo algorithm for 
lattice trees J. Phys. A: Math. Gen. 25 (1992) 303-333. 

[26] E. Orlandini, F. Seno, A.L. Stella and M.C. Tesi, Collapse from linear to branched 
polymer behavior, Phys. Rev. Lett. 68 (1992) 488-491, see also L. Turban, Gener- 
alised self-avoiding walk, J. Phys. A: Math. Gen. 16 (1983) L643-646. 



36 



